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We discuss the phase structure of a lattice Higgs- Yukawa system in the variational 
mean field approximation with contributions of fermionic determinant being calculated 
in a ladder approximation. In particular, we demonstrate that in this approximation the 
ferrimagnetic phase in the Z2 model with naive fermions can appear as an artifact of a finite 
lattice and that the phase diagram for this model on infinite lattice changes qualitatively 
at space-time dimension D = 4 compared with those at D > 4. 
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1. Introduction 

Although mean field method for lattice systems including fermions loses consider- 
ably its simplicity and requires further approximations it is still useful to get some idea 
of the phase structure of the systems and to orientate Monte Carlo simulations towards 
investigating the most interesting points. In this paper we make an improvement in the 
approximations within the variational mean field approximation for Z2 Higgs- Yukawa sys- 
tems by summing up a ladder type contributions to fermionic determinant, including those 
of the next order in inverse space-time dimension 1/D. This enables us to observe two new 
points. As the first one we demonstrate that within our approximation the ferrimagnetic 
phase in the simplest Higgs- Yukawa model with naive fermions can arise as a finite lattice 
artifact. The second point is that the value D = 4 turns out in a sense to be critical, as 
the domain of paramagnetic phase just at D = 4 becomes disconnected, being connected 



The paper is organized as follows. The system under consideration is defined in Sect. 
2. In Sect. 3 we describe the method and approximations. Results are discussed in Sect. 4. 

2. The model 

The system is defined on a hyper cubic D-dimensional (D is even) lattice A with sites 
numbered by n = (ni, ...,n_o), — N/2 + 1 < < N/2 (N is even) and with lattice spacing 
a = 1; /* is the unit vector along the lattice link in the positive /(/-direction. Dynamical 
variables of the model are the fermion 2 r> / 2 -component fields ip n , ip n , and scalar field 
4>n £ Z2 (i.e. 4> n = ±1). We imply antiperiodic boundary conditions for the fermion and 
periodic for the scalar fields. 

The model is defined by functional integral 



at D > 4. 
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with the action 
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K G (— oo, oo) is the hoping parameter, y > is the Yukawa coupling; we use the Hermitean 
7-matrices: [7^,7^]+ = 28^; L^(p) = smp^, p^ = (2-n/N)(k^ - 1/2), -N/2 + 1 < k^ < 
N/2, so that p M G (— tt/2, 7r/2). Operator satisfies 

0mn = 0nm- (2-4) 

In the limit of N — > oo the sum iV _r) defines the integral J^^ 2 d D p/ (2n) D . 
The action (|2j|) is invariant under Z2 global chiral transformations 

0n^-0n, ^^(-Pl+Pr)^, ? n ~> ("Pi? + P L ) , (2.5) 

where = (1 ± jr>+i)/2 are chiral projecting operators. 



3. The method and approximations 

To analyze the phase structure of the model we use the variational mean field ap- 
proximation |]J (see also Pj) which becomes applicable to ( |2.1j ) after integrating out the 
fermions 

-W[J] ^ 2 «En, M ^ 

+ In det [0 + y4>] + J2 n J n (p n 

Z [ J \ = e = 2^ e • t 3 - 1 ) 

4>n&Z 2 

Then for free energy of the system F = W[0] the method yields the inequality 

F < F MF = inf [- J^M/in) - Ku'(h n )) - (2k £ n n+A + In det [0 + j/0]) h ], (3.2) 

n n,fi 

where hr, is a mean field, and 



u(h n ) = In e n ^ n =ln2cosh/i n , 
(O[0]), = e" E ^ W £ O[0]e E ^ r 



(3-3) 



So, we can get some idea of the system, studying Fmf, that is much simpler than that for 
F. From ( |3.3j ) it immediately follows that 



t ) = u'{h n ) = tanh/i n , 

(3.4) 

n4>n)h = u'(h m )u'(h n ) + 5 mn u"(h m ), etc., 
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and therefore the main problem is a calculation of the expectation value of the fermionic 
determinant 

(lndet[0 + y(/>]) h 

^ (-1) 



= lndet0]-^^-^y" £ ^fe^l-CD^^-^). 



n 



(3.5) 



f — ll n 1 

= 2°/ 2 iV D ln 2/ -^^-^- £ tr^.Aia-^nu)^^-^)^ 

n=l 

where tr stands for the trace over spinorial indices; in the first term of the second equation 
the relation 0? 1 = 1 has been taken into account. 

Following the usual way we consider Fmf for two translation invariant ansatzes for 



K M = h, 

,AF _ , u , -( ^\L i^ n » 



(3.6) 



K* = e n h, e n = (-I)' 

which in fact are the order parameters distinguishing the ferromagnetic (FM: h FM ^ 0, 
h AF = 0), antiferromagnetic (AF: h FM = 0, h^ F ^ 0), paramagnetic (PM: both are 
zero), and ferrimagnetic (FI: both are nonzero) phases in the system. Then the mean field 
equations are reduced to 

9 jpFM, AF _ n ,o 7 \ 

~Qy l t MF - u ' O'J 

where F^^ ,AF is the functional of the right-hand side of Eq.(3.2) on ansatzes (3.6). Fur- 
ther simplification comes from the observation (see, for example 0), that as the value 
h = is always a solution of Eq. (3.7), and, therefore, second order phase transition lines 
are determined by equations 

° pFM.AF I _ n /o o\ 

gj^^MF \h=0- U, (3.S) 

to find them it is sufficient to know (lndet [0 + y<p])h to terms of order of h 2 . 

If the problem could be solved exactly both of two representations (3.5) of the fermionic 
determinant would yield the same answer. But correlations of </>'s at coinciding arguments 
(Eq.(3.4)) make the problem unsolvable exactly, as the contributions of order of h 2 to (3.5) 
come from terms of any orders of u" , as well as from those of order of u' 2 . These contribu- 
tions shown schematically in Fig. 1. Therefore, we are forced to use some approximations, 
and, particularly, to use two representations of (3.5) separately for "weak" and "strong" 
coupling regimes of y, though the exact meaning of this can only be clear a posteriori. 
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Our approximation involves summing up all diagrams of Fig.l (a) (proper ladder di- 
agrams) and (b) (crossed ladder diagrams), so we may call it as a ladder approximation. 
Using property (2.4) of the Dirac operator we find that the contributions to F™' AF from 
the fermionic determinant, AF^^f' AF , have the same functional form for both representa- 
tions (3.5) and in our approximation read as follows 

a rpFM AT D D/2-l, C 2 u' 2 G(0) D \- (?u"G{q) 

AF MF=N2 ( 1 + c 2 u , fG{Q) + N L, 1 + c 2 u , G(q) ) 

a -pAM at DqD/2-u c 2 u' 2 G(n) D \- c 2 u"G(q) 

AF MF =N2? ( 1 + c2u „ G{iT) +N ^ 1 + c2u „ G{q) ), 

where g M = (2ir/N)l fJi , —N/2 + 1 < < N/2 (so that g M G (— 7r, 7r]), while coupling c and 
form of function G depend on the representation. So, for weak coupling regime we have 
c = y and 

and for the strong coupling they are c = y~ x and 

G S (q) = J^dY1 L (P) L (P + 9) = E sin ^ sin ^ + ( 3 - n ) 

The first terms in (3.9) come from the diagrams of Fig. 1(a), while the second from those 
of Fig.l(ft). 

Then, from Eq. (3.8) and the above formulae it follows that critical lines in the system 
in our approximation are determined by the expressions 



K F(W) = J_M _ qD/2 2/ G W (0) _ D G W (q) 

AD [ g{ l + y 2 G w (0) Y {l + y 2 G w {q)) 2)h 

G W (n) AT _ D ^ G w (q) 



af(w) Lh _ 2 D / 2 v 2 ( ° ^ N~ D V — 

Kcr ~ 4D [i 1 V[ l + y 2 G^M ^(1 



l + y 2 G^(n) ^ (l + y 2 G^(q)) 2 



(3.12) 



4D L V + G s (0 ^ (y 2 + G s (q)¥ 



af(S) Lh _ 2 D/2 ( G ^ - N~ D V V ° ^ )] 

4D [ { y 2 + G s {7r) ^ (y 2 + G s (q)) 2)l ~ 

We now should make some comments. 
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(i) The contributions to (3.12) which are proportional to G(0) and G(tt) are general- 
ization of "double chain" contributions of Ref. @, as the diagrams of Fig. 1(a) are the gen- 
eralization of the double chains to any configurations of the same topology. They coincide 
only for G s because of strict locality of the Dirac operator, but not for G w . More impor- 
tant difference comes from the second terms corresponding to the diagrams of Fig. 1(6) (the 
latter correspond to the generalization of the double chains with coinciding ends), which 
have not been taken into account in previous calculations (see also 0). From the well 
known symmetry of the model under the transformations: (if), if)) n — > exp(ze n 7r/4) (■*/>, if)) n , 
4>n — » ^n4>ni k — ► —k, y — > — iy, it follows that G(tt) = — G(0), and also, that the contribu- 
tions of the new terms are of even power in y ±2 beginning from y ±4 . 

(ii) These terms can become dominating when y 2 is close to the values 1/G W (0) 
or G s (0) which are singular points of the expressions under the sum, even though in 
weak coupling regime they are of 0(D~ l ) compared with the first ones. Thereby these 
terms determine domains of the "weak" and "strong" coupling regimes also for n^ r . They 
are domains of analyticity of functions n^.(y) and K^ r (y), that is y 2 < 1/G W (0) and 
y 2 > G s (0), respectively, coinciding for and K^ r F . 

(iii) We have no strict arguments why other diagrams which we did not take into 
account could be neglected compared with the ladder ones. In particular, in strong coupling 
regime they can give contributions to n cr of the same order in 1/D as the latter. But 
because those diagrams come into play in higher orders in y , at least from the order 
of y ±6 , the assumption that their contributions are suppressed and less singular looks 
plausible. 

Finally, it worth noting that the formulae (3.12) are applicable to any lattice fermion 
actions, including non-local ones, whose Dirac operators satisfy property (2.4) Q. 

4. Results and discussion 

Let us now compare the phase diagrams determined by the expressions (3.12) for 
D = 4 for finite iV and for the limiting case of iV — > oo . The new terms are always 
negative and therefore increase the contributions of the first terms for and decrease 
them for k^JF . The question is of how much. 

For JV = 4 we have G w (0) = 0.5, G s (0) = 2, so that the domain of inapplicability 
of our formulae shrinks to the point y = 2 1 / 2 , and the phase diagram is shown in Fig. 
2. The curves K>f r (y) and K^F(y) intersect each other forming narrow domain with the 
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ferrimagnetic phase around y = 2 1 / 2 , which is spreaded from — oo to oo in k. It is natural to 
assume, that contributions of other diagrams (Fig. 1(c)) smooth the negative contribution 
of those of Fig. 1(b), so that the PM-AF phase transition line in Fig. 2(6) becomes 
continuous. Then, as a result we would have a familiar picture, typical for SU(2) models 
(see, for example, ||), with FI phase lying below this line. 

In the limit of N — > oo we have G w (0) ~ 0.62, G s (0) = 2, but the picture is changed 
qualitatively. The phase diagram is shown in Fig. 3. The curves do not touch each other 
even at k — > — oo and FI phase does not appear. 

To clear up why this happens let us consider behaviour of functions of y 2 determined 
by the sums in (3.12) near the points 1/G W (0) and G s (0). Let us define positive 5 = 
y 2 — G s (0) or 1/G W (0) — y 2 . Then a simple analysis shows that at a finite N and a small 
8's these functions is of order of 0(D k N~ D 5~ 2 ), so the intersections of the curves K^ r (y) 
and K^ r F (y) always occur at the points 5 = 0(D l N~ D ), -k = 0(D m 2 D / 2 N D ), where k, 
/, m are some (negative or non- negative) powers. But at N —>■ oo, when the sums go over 
to integrals, this functions become of order of In S for D = 4, and even of O(S ) at D > 4. 
This means that at D > 4 we can continue lines Kcr\y) and Kcr (y) until they intersect 
each other, so that the phase diagram in this case looks like in Fig. 4, that reproduces the 
result of ref. . 

Thus, this example demonstrates importance of summing up contributions to 
fermionic determinant including those of the next order in 1/D for D — 4 systems. An- 
other point is that even though we cannot definitely conclude whether the FI phase in this 
example is an artifact only of a finite lattice or also of the mean field approximation, this 
gives one one more caution in what concerns finite lattice effects. 
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Fig. 1. Diagrams contributed to AF to order h 2 . Solid lines denote or 1 , each 
solid circle stands for it', dashed line for u" . 
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Fig. 2. (a) Phase diagram of the model at D = 4, N = 4. Intersections of FM-PM 
phase transition line (solid) with PM-AF phase transition line (gray) form FI phase in the 
narrow region around the point y = 2 1 / 2 shown in (6) in more detail. 
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Fig. 4. Qualitative picture of the phase diagram of the model at D > 4, N — > oo. 
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